LINEAR REPRESENTATIONS AND ISOSPECTRALITY WITH 
BOUNDARY CONDITIONS 



ORI PARZANCHEVSKli AND RAM BAND^ 

Abstract. We present a method for constructing families of isospectral sys- 
tems, using linear representations of finite groups. We focus on quantum 
graphs, for which we give a complete treatment. However, the method pre- 
sented can be applied to other systems such as manifolds and two-dimensional 
drums. This is demonstrated by reproducing some known isospectral drums, 
and new examples are obtained as well. In particular, Sunada's method [Ij is 
a special case of the one presented. 



1. Introduction 

"Can one hear the shape of a drum?" - This question was posed by Marc Kac in 
1966 [2]. In other words, is it possible to determine the shape of a planar Euclidean 
domain from the spectrum of the Laplace operator on it? This question gave rise 
to fertile research, investigating it from various aspects. Two main approaches 
were, on the one hand, attempts to deal with the inverse question of reconstructing 
the shape from the spectrum, and on the other hand, trying to find systems whose 
shapes are different, yet have the same spectrum. Such examples are called isospec- 
tral. Although Kac's original question regarded two dimensional planar drums, the 
research on isospectrality expanded quickly to other types of systems. We will not 
go into detail, but refer the interested reader to [l]-[9] for a broader view of the 
field. However, we will mention here two milestones in the field of isospectrality. A 
theorem by Sunada gave an important machinery for the construction of isospec- 
tral Riemannian manifolds [Ij. Later, this method was used by Gordon, Webb and 
Wolpert to construct the first pair of isospectral planar Euclidean domains [21 [4] 
thus negatively answering Kac's original question. 

This paper starts with a presentation of the basic theory of quantum graphs and 
existing results on quantum graph isospectrality. We then present the algebraic 
part of our theory and its main theorem. This is followed by a section which 
explains the construction of the so called quotient graphs that lie in the heart 
of the theory. After the theory is fully presented, we apply it to obtain various 
examples of isospectral quantum graphs. We then demonstrate how to apply the 
method to other systems, explaining some known results, as well as obtaining new 
ones. In particular we discuss the relation to Sunada's method. We conclude by 
pointing out key elements of the theory that are to be investigated further and by 
presenting open questions. 
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2. Quantum graphs 

A graph T consists of a finite set of vertices V = {vi} and a finite set E = {cj} of 
edges connecting the vertices. We assume that there are no parahel edges (different 
edges with the same endpoints) or loops (edges connecting a vertex to itself), but 
we shall see that this infiicts only a small loss of generality. We denote by Ey the 
set of all edges incident to the vertex v. The degree (valency) of the vertex v is 
dy = \Ey\. T becomes a metric graph if each edge e e E is assigned a finite length 
le > 0. It is then possible to identify the edge e with a finite segment [0, le] of the 
real line, having the natural coordinate Xe along it. A function on the graph is a 

vector / = (/|ei ' • • • ' /|e|£;|) functions : [0,/ej] C on the edges. We shall 
usually consider smooth functions on the graph, meaning that f\^^ C°°([0, le]) for 
all e e E. Notice that in general it is not required that for v G V and e, G Ey 
the functions /|^ and f\^, agree on v. 

To obtain a quantum graph, we consider a differential operator on the graph, 
by default the negative Laplacian: —A/ = ^— /''|^^ , . . . , — /'^^^ J . In addition, 

we require the functions on the graph to obey certain boundary conditions stated 
a priori; for each vertex v G V, we consider homogeneous boundary conditions 
which involve the values and derivatives of the function at the vertex, of the form 
' f\^ -\- By • f'\^ = 0. Here Ay and By are dy x dy complex matrices, / 



IS 



the vector (^/L^/^) ' 



and 



/|e^^ i^)^ of the values of / on the edges in Ey at v, 

' "^'le^d '^^^) i^ vector of outgoing derivatives of 
/ taken at the vertex. To sum up, a quantum graph is a metric graph equipped 
with a differential operator and homogeneous differential boundary conditions at 
the vertices. Notice that before stating the boundary conditions, the graph is 
merely a collection of independent edges with functions defined separately on each 
edge. The connectivity of the graph is manifested through the boundary condi- 
tions. We denote by C°°(r) the space of (smooth complex) functions on the graph 
which satisfy the boundary conditions at the vertices, and by H (F) the subspace of 
C°°(r) spanned by eigenfunctions of the Laplacian. The reader interested in more 
information about quantum graphs is referred to the reviews fTO \ \TT\ \T2\. 

A standard choice of boundary conditions which we adopt is the so called Neu- 
mann boundary conditiorQ: 



/ agrees on the vertices: e V Me^e' ^ Ey \ f {v) = f {v). 



• The sum of outgoing derivatives at each vertex is zero: G F : X] f 
The Neumann boundary condition can thus be represented by the matrices 



{v) = 0. 



Ay 



/I 



-1 



/O 



-1 

0/ 



By 









1/ 



For a vertex of degree one the Neumann condition is expressed by the matrices 
Ay = (O), By = (l), and means that the derivative of the function is zero at the 



■"^This condition is also widely encountered under the name of Kirchhoff condition. 
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leaf V. Another natural boundary condition for leaves is the Dirichlet boundary 
condition: Ay = (l), By = (o), which means that the function vanishes at the 
vertex. 

Neumann vertices of degree two deserve a special attention. They can be thought 
of as inner points along a single edge - the concatenation of the two edges incident 
to the vertex - and we would like to be able to add or remove such inner points, for 
reasons which will become clear later on. At such points, however, a function on the 
edge is only required (by the Neumann condition) to be continuously differentiable 
(C^), rather then smooth (C ); therefore, adding a Neumann vertex of degree 
two at an inner point of an edge augments the space of allowed functions (by ones 
such as 1 3y I ' X ). The question of C versus C is inherent to the modeling of 
one dimensional manifolds as quantum graphs. For example, in order to regard the 
circle as a quantum graph, we must place at least one vertex along it, and at this 
vertex functions on the resulting graph may have a non-differentiable derivative. 

The good news is that adding or removing Neumann vertices of degree two 
does not change the spectral properties of the graph in question. For sums of 
eigenfunctions of the Laplacian, being and piecewise C is equivalent to being 
altogether, so that if the graph F' is obtained from F by adding or removing 
such points, we have H {T') = H (F). With this observation in mind, we will allow 
ourselves to make manipulations of degree two Neumann vertices, with no essential 
loss of generality (at least from the spectral viewpoint). For example, loops and 
parallel edges can be eliminated by the introduction of such "dummy" vertices, so 
that as mentioned, we shall assume that we are dealing with graphs with no such 
nuisances. 

If for every v eV the dy x 2dy matrix {Ay \ By) is of full rank, we shall say that 
the quantum graph is exact Non-exact quantum graphs are not very interesting 
from the spectral point of view, as their spectrum is all of C. On the other hand, 
we shall later be led to consider the opposite phenomena, i.e., vertices at which 
there are "too many" boundary conditions. In this case we shall admit Ay and By 
to be of size m x dy, possibly with m > dy, and we shall call the corresponding 
graphs generalized quantum graphs. From the spectral perspective these are much 
more interesting than non-exact quantum graphs. Consider for example a Y-shaped 
graph, with a Neumann condition at the center, Dirichlet conditions at two of the 
leaves, and the condition Ay = (q), By = ( ? ) ^it the third; its spectrum is nonempty 
if and only if the lengths of the two edges with Dirichlet leaves are commensurable. 

There is a natural inner product on C°°(F), given by {f,g) = /o"^ f\ ' ^^e- 

eeE ^ ^ 

Kostrykin and Schrader [13j provide necessary and sufficient conditions for the 
Laplacian to be self-adjoint with respect to this product. These conditions can be 
stated in a number of equivalent forms (see [IS]). We give two of them: 

(1) F is exact, and Ay • Bl is self- adjoint for every v eV. 

(2) For every v e V there exist a unitary matrix U such that {Ay \ By) is 
row-e quivalent to {i{U - I) \ U ^ I) 0. 

In particular, Neumann and Dirichlet boundary conditions satisfy these require- 
ments. 

There are several known results concerning isospectrality of quantum graphs. 
Gutkin and Smilansky [15] show that under certain conditions a quantum graph 



^(2) follows from (1) by taking U = 2 {-A + iB) M + /. The other direction is trivial. 
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can be heard, meaning that it can be recovered from the spectrum of its Laplacian. 
On the other hand, constructions of isospectral graphs were also estabhshed, by 
various means: by a trace formula for the heat kernel [16j, by turning isospectral 
discrete graphs into equilateral quantum graphs [17J, and weighted discrete graphs 
into non-equilateral ones [H] ; in [T5t [T9] a wealth of examples is given by an analogy 
to the isospectral drums obtained by Buser et al. [5], and in [20j is presented an 
example, whose generalization has led to the theory presented in this paper. 

3. Algebra 

For a quantum graph F we can regard H (F) as a C [x]-module, where x acts as 
the (negative) Laplacian. For every A G C we denote by <l>r (A) the submodule 

$r (A) = Ann„(r) (x - A) = {/ e W (r) | - A/ = A/} , 

which as a vector space is merely the A-eigenspace of the Laplacian. The spectrum 
of F is the function 

(Jr : A dime (A) , 
which assigns to each eigenvalue its multiplicit}|l. Two quantum graphs F and F' 
are said to be isospectral if their spectra coincide, that is cjy = crp', and as noted 
in [21j, this can follow from the stronger assumption that H{T) and (F') are 
isomorphic as C [x]-modules, which means that their Laplacians are conjugate. 

A symmetry of a quantum graph is an invertible graph map that preserves both 
the lengths of edges and the boundary conditions at the vertices. The group of 
all such symmetries is denoted AutF. A left action of a group G on a quantum 
graph F is equivalent to a group homomorphism G AutF. Such action induces 
a left action of G on H{T) (by (gf) {x) = f [g~^x) - the inversion accounts for 
the contravariantness of H). This gives H (F) a CG [x]-module structure, since the 
Laplacian commutes with all symmetries. The eigenspaces $r (A) = Ann {x — A) 
are again submodules, and in particular they are CG-modules, that is, complex 
representations of G. Assuming that G is finite, with irreducible complex represen- 
tations /Si, . . . , 5'r, we can decompose each eigenspace to its isotypic components: 

r 

(3.1) $r(A) = 0*^*(A) , 

i=l 

where (A) = S'i . . . S'i as CG-modules. 

We start by counting separately, for each irreducible representation S of G, only 
the A-eigenfunctions which reside in W- This means that we are restricting our 
attention to functions which under the action of CG span a space that is isomorphic, 
as a representation of G, to S. However, since dimS* always divides dim^p (A), we 
can already normalize by it. We thus define the spectrum of S as 

(3.2) cr^ : A dime <3E>F(A)/dimc S . 

By the orthogonality relations of irreducible characters, we can rewrite this as 
cTp (A) = (xs-, X^r(A))c" expanding this linearly, we define the spectrum of 
for every representation R of G, to be 

(3.3) cr?:Ah^(xi?,X<i>r(A))G • 



^In effect we have ot : C ^ {0..2 as the eigenvalue of a Laplacian eigenfunction, together 
with the values {/|g (0) , f'\^ ^^^}eeE^ determine the function. 
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cTp (A) has an algebraic significance: it refiects the size of the S'-isotypic part of 
$r (A). Looking for a parahel algebraic interpretation of (A), we find that 

(A) = dime HomcG (A)) . 

Quite generally, if A ^ 5 is a ring extension and M and N are modules over A 
and B respectively, then for Cb {A), the centralizer of A in 50, HomA (M, N) has 
a natural Cb (A)-module structure (by {bf) (m) = b ■ f (m) for every b e Cb (A)). 
For our purposes, since C [x] lies in the centralizer of CO in CG [x], we obtain that 
HomcG {R^ ^ (r)) has a C [x]-module structure: 

X • / : r ^ - A (/ (r)) (/ G Homcc {R. H (F))) . 

This allows us to make the following definition. 

Definition 1. A ^/n-graph is any quantum graph T' such that there is a C[x]- 
module isomorphism 

(3.4) n{T')=RomcG{R.n{T)) . 

We note, in particular, that for such F', there is an isomorphism 
$r'(A) = 

Ann'^(r/) {x - X) ^ AnnHomcG(i?,'H(r)) (^ - A) 
= HomcG (^,^r (A)) 

which by taking dimensions translates to equality of spectra: 

(3.5) err' = cr^ . 

Since is not a spectrum in the classical sense, we cannot really call this 
isospectrality. However we do have from this that all r/i^-graphs are isospectral 
to one another, and we will use this to speak non-rigorously about "the spectrum 
of ^/r^\ (Jr/ji = cTp . The following proposition exhibits another manifestation of 
isospectrality. 

Proposition 2. All ^/cG-graphs are isospectral to F. 

Proof. By (|3.ip . (|3.2p . and linearity, the classical spectrum ar coincides with the 
spectrum of the regular representation of G: 

r 

(3.6) (7^^ = dim Si • crf ^ = ar . 

i=l 

This can also be deduced from the fact that for every i?-module M there is an 
isomorphism Roirr {R, M) = M, so that we have HomcG (CG, (A)) = (A) 
for every eigenvalue A. □ 

We can say even more: 

Theorem 3. Let T be a quantum graph equipped with an action of G, H a subgroup 
of G, and R a representation of H. Then ^/r is isospectral to r/indgi?. 



^That is, Cb {A) = {h ^ B\ya ^ A : ah = ha}. 
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Proof. This follows at once from the Frobenius Reciprocity Theorem, which states 

that there is an isomorphism Home// (i?, (F)) = HomcG (ind^R^H {T)^ . It is 

straightforward to verify that this is an isomorphism of C [x] -modules. Note that 
from the formal point of view, we have actually shown that ^/r and ^/ind%R are 
identical (as classes of quantum graphs). □ 

Remark. This gives yet another explanation for the equality of the classical spec- 
trum with that of the regular representation (proposition [2]) : for H — and 
\h its trivial representation, it is clear by the isotypic component perspective that 
{n (F))^^ = n (F), so that jSSl) follows from Indglj:/ ^ CG. 

Corollary 4. IjG acts on F and i^i, H2 are subgroups of G with corresponding 
representations R2, such that Ind^^i^i = Ind^2^2; then ^/Ri and ^/r2 are 
is spectral. 

Remark. This corollary is in fact equivalent to the theorem, which follows by taking 
H2 = G, R2 = Ind^^i^i. It is presented for being of practical usefulness (it allows 
one to work with representations of lower dimension, as can be seen in section [5]), 
but also since it indicates the bridge connecting our method with the classical one 
of Sunada. In section [631 we shall cross it. 

The sharpest observations in this section would be mere algebraic tautologies, 
unless we can show that r/i^-graphs do exist. The next section is devoted to this 
purpose. 

4. Building t/^-graphs 

In this section we prove the existence of the quotient graphs ^/r. This is done 
by describing an explicit construction of '^/r, given a graph F, a representation 
R of some group G acting on the graph, and various choices of bases for this 
representation. As the lengthy technical details of the construction might encloud 
the essence of the method, the reader may prefer to go over sect ion [5]first, and obtain 
an intuition for the construction of the quotient graph from the examples presented 
there. More intuition for the construction can be gained from the examples in |14| . 

We summarize the main conclusions of this section in the following theorem: 

Theorem 5. For any representation R of a finite group G, which acts upon a 
quantum graph T, there exists a generalized ^/r quantum graph. Furthermore, if 
T^s Laplacian is self-adjoint, then there exists a proper ^/r quantum graph, and it 
is exact. 

4.1. Intuition. A motivation for the construction of our quotient graphs is given 
by thinking about it as an "encoding scheme" 0. An element / in HomcG (^, ^ (r)) 
can be thought of as a family of functions on F, parametrized by R. To emphasize 
this view, we shall write fr instead of / (r) (where r e R). Our goal is to build 
a new graph, each of whose complex functions encodes exactly one such family. 
The desired map ^ : HomcG (^, (r)) HC^/r) (see definition [T]) is in fact 
this encoding. An encoding scheme should always be injective (in order to allow 
decoding), but we have also required ^ to be surjective: this can be translated 



^In section [6] we show that the same construction and motivation can be appHed analogously 
to other geometric systems. 
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to the idea that the encoding must be "as efficient as possible" |j - that ^/r is 
to be a "minimal" graph allowing such an encoding, since it admits no Laplacian 
eigenfunctions apart from the ones used by the scheme. 

First, we reduce the infinite family / to a finite one by choosing a basis B — 

{hj} .^-^ for and restricting our attention to I fb- > . From these "basis func- 

^ ^^^^ 
tions" we can reconstruct /, since the CG-linearity of / implies in particular C- 

linearity (i.e., f^ajbj = ^(^jfbj)- As a first encoding attempt we could take a 

graph with d times each edge in F, and let the j*^ copy of the edge e carry the j*^ 



fbj 



However, this en- 



basis function restricted to e. That is, define (^^/^ 
coding is not efficient enough, since we have used only C-linearity. For each ^ G G, 

determines < fr > , specifically by 



eR 



CG-linearity implies that < fJ \ 

I I e J r 

(4.1) fr\ ^ (g-'fr 



(the inversion occurs since G acts on (F) by g ■ f = f ^ g~^)- Thus, it suffices to 
encode the basis functions on only one edge from each G-orbit of F's edges. 

It turns out that if the action of G on £^ is free, then apart from determining 
the appropriate boundary conditions at the vertices we are done: for {e*}, a choice 



of representatives for ^/g, setting (^^/^ 



fbj 



(where 1 < j < <i) is indeed a 



"good" encoding (i.e., once the boundary conditions are correctly stated, ^ is an 
isomorphism.) 

If, however, some edge e — {y^v'} has a non-trivial stabilizer Gg = G^ H G^/, 
then greater efficiency can (and therefore must) be achieved. For example, assume 
that dimi? = 1 and that for some ^ G Ge we have g ^ kerp^^, where pR is the 
structure homomorphism G GLi (C). We then have 



/r 

which implies that fj 



pRi9)-r 



^ PR (g) • fr 



for all r, and as a result, the edge e need not have any 



representative in the quotient. We can "decode fr\ from thin air", since we know in 

I e 

advance that it can only be the zero function. The generalization of this observation 
is that for each edge e, the information in < > is encapsulated in R^"" 0: if 

le J reR 

r belongs to a nontrivial component of Res^^i?, then fr = 0. Therefore, we need 

only di = dimi^^ei copies of each representative in the quotieniH. This further 
"compression" slightly complicates the determination of the boundary conditions. 
When G acted freely on the edges, we had d functions, fi)., each satisfying the 
boundary conditions at the vertices of F, and we could have translated this quite 
easily to boundary conditions on the quotient. Now, however, for each edge we 



^In a suitable sense, since better encoding may exist, but we want the encoding to be by 
another quantum graph, in a manner which intertwines the corresponding Laplacians. 

^R^ is the trivial component of Res^R, i.e. R^ = {r e R\yh e H : hr = r}. 

^However, we shall later find it convenient to think about d = dim R copies, where the di+l . . .d 
copies are "dead", meaning that whenever a function on them appears in a formula it is to be 
understood as zero. 
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need to encode a "function basis" j/b* | where {6* } is a basis for R^^K Since 

for different e*'s the spaces R^^^ need not even overlap, we now have only function- 
chunks, indexed by different i?-elements for each edge, and no function on the whole 
of r to extract boundary conditions from. Fortunately, algebra is generous and this 
complication turns out to be solvable. 



4.2. Method. We now present the actual construction procedure. Assume we have 
a representation of a group G acting on the quantum graph T = V), and we 

have chosen representatives {e*} .^-^ for the orbits ^/g, and likewise {vk}^^i for ^/g. 

We have also chosen an ordered basis B = (hj)^.^-^ for i?, and for each i e {1../} 

another ordered basis for R, = such that is a basis for R^^^ 

and each b] with j > di = dim R^^' lies in a nontrivial component of Res^ . R. 

The quotient graph ^/r obtained from these choices is defined to have {vk}j^=i 

as its set of vertices, and {e^}*"^ ^ for edges, where each ej- is of length /gi. If 

connects gvk to g'vk' in F, then, for all j, connects Vk to Vk' in ^/r. We shall 
assume, by adding "dummy" vertices if needed, that G does not carry any vertex 
in V to one of its neighbors. This serves three purposes: 

(1) It means that ^/r has no loops; i.e., that k ^ k' m. the notation above. 



This allows us to speak of / 



(v), the value of at without confusion 



regarding which end of e* is meant. 

(2) It assures that an edge is not transformed onto itself in the opposite direc- 
tion, in which case we would have had to take only half of the edge as a 
representative for its orbit. 

(3) It assures that the fixed points of each g ^ G are either entire edges, or 
vertices. 

Note that in order that G still act on the graph, the dummy vertices are to be 
added in accordance with its action, i.e., if a vertex is placed at x G (0, le) along e, 
one should also be placed at x along ge, for every g e G. 
We can now define ^ on HomcG {R, H{T)): 



de_f ~ 



and it is clear that ^ does intertwine the Laplacians. We would like to determine 
vertex conditions on ^/r that will ensure that ^ is into and onto H {^/r). This will 
require nothing more than linear algebra, and we start by rephrasing (|4.ip basis- 
wise. We make the following convention: an expression in bold is to be understood 
as a row vector of length d, where the ^-symbol indicates the place of the index; 

e.g., stands for (^fbi\^, •.• , /bJe)- 

Consider fr , an arbitrary function in the family / evaluated on an arbitrary 
edge. Write r G as 6 • a, where b = . . . , 6^^) G Mixd {R) and a G Mdxi (C). 
r = b ■ cx implies gr = b\ ■ [pR {g)]%i ■ a, and therefore, by (|4.ip we have fr = 



LINEAR REPRESENTATIONS AND ISOSPECTRALITY WITH BOUNDARY CONDITIONS 9 



Linearity now implies 



where pR is understood, f^. 



J • ' • ^ fbi 



/j, and 

,0,...,0 I , 



since for j > di we have seen that /^i =0, and we therefore did not include the 

corresponding edge in ^/r (it is "dead" - see footnote [8]) . We now see that for 
/ G (V^) inverse of ^ must be given by: 



hex. 



(again stands for [pr {9~^)] ^i), so we need to establish that the r.h.s 

is independent in the choice of g. We observe that if g and g' are two possible 
choices then g~^g' G G^m ^nd by the construction of we have [9~^9']^i = 





^ 





* 



. As we have agreed that / 



/ 



, 0, . . . , , we have 



and thus 



/ 



B 



establishing that ^ ^ is well defined. 

We can now determine matrices Ay^ and By^ for the vertex Vk from the matrices 
Ay^^ By^ of the vertex Vk- Assume that the edges entering Vk are ^le^S . . . ^gn^^"^ 
(where n = dy^)^ so that a function / on F satisfies the vertex conditions at Vk 
when 



where we recall from section [2] that 



0, 



f 



(/' . fe) 



r 



/ G HomcG {R^H (r)) means that fr satisfies the conditions at Vk for all r e R, 
which happens iff the basis functions | satisfy them. Thus, if we define the 
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X d matrix 



ft 



( hi 



fb. 



and analogously 
(4.2) 



then we need only check that 



' fb 



-^Vk ' fb 



In addition, we note that if the boundary conditions are met by / at Vkj then 
they are also met at any vertex in the orbit G - Vkj since G is assumed to preserve 
boundary conditions. 

For a n X m matrix X = {{xij)) we define its row- wise- vectorization to be the 
nm X 1 matrix 



def 



(xxx 7^12? • • • 1 Xirni ^21 ? ? "^nm J 

\ (yXnl , . . . , ^nm) 

Vectorization behaves quite nicely under multiplication. Specifically, ty {A - B - C) 
[a (g) C^) ■ TV B, which allows us to write (|4.2p as 

(4.3) {Ay, Id) • rv _ + {B^, Id) • rv _ = On, 



dxl 



Recalling that ft^ 



^f 



' [Si we have 



\ 



( 



9i 



and likewise for rv/^ 



But now, the last vector is almost ^/ 



the vector 



of values of (^^/^ ^/c' Only two changes need to be made: first, if the edges 
entering Vk are . . . , e^^ , e^^, . . . , e^"^ , then by definition {/ii}^x = {i^ij^^^ 
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sets; however, the /i^ are distinct, whereas in general, repetitions can occur among 
the Vi (i.e., two edges in Ey^ might belong to the same G-orbit). Second, as in 
all our expressions there might be "dead" edges, e^^ with j > d^. , which do not 
really appear in the quotient graph (note, however, that neither of the problems can 
occur when the action of G is free). We shall deal with these two inconveniences 

1 Ui = jlj 



and then take 



at once: we define the n x m matrix (6^) - • = 

I otherwise 

O to be the nd x matrix obtained by removing from (0' Id) the columns 
{{i — 1) ■ d -\- j} i<i<m ; these are the columns which would have been multiplied 

<j<d 

by a "dead" edge in f ^/L^i ••• ^/L^i ^f\,^^2 ^/L^m ) . We now have 

\ 1 d 1 d / 



6 • ( ^i^^'^) 



) =6-^/ 



and we can thus define 

(4.4) 
(4.5) 



where = diag {^g^ ? • • • ? [^n ^l^^r.) ' finally rewrite (|4.3p as 



= 0, 



These vertex conditions on ^ f at Vk are equivalent to fr satisfying the vertex 
conditions at Vk for all r e and therefore also on the entire orbit G - Vk- If we 
repeat this process for each /c = 1 . . . K, we indeed obtain boundary conditions on 
r/i? which are satisfied by ^/ exactly when / G HomcG (^, H (r)). 

If the action of G is free, then is just a permutation matrix (we can even order 
Eyj^ SO that 9 = /) , but in the general case 6 might be non-square (explicitly, it is 
of size ndxdy^^ where dy^ = Yl^i ^fj'i — — ^^)- When this occurs, the matrices 
Ayj^ and By^^ we have obtained are not square matrices, and we therefore obtain a 
quotient which is only a generalized quantum graph. Nevertheless, as the matrices 

0, we can perform elementary row operations on the nd x 2dy^ matrix [Ay^ \ By^) 
without changing the boundary conditions at Vk^ and thus perhaps reduce the 
number of rows of {Ay^ \ By^ ) . In the case that rank {Ay^ \ By^ ) < dy^^ we can reduce 
the matrices Ay^ and By^ to squares ones, and if this holds for all k then we actually 
have a proper quantum graph. If it further happens that rank {Ay^ \ By^) = dy^ for 
all /c, then the quotient graph is also exact. We now show sufficient conditions for 
this to happen. 

Proposition 6. // there exist G and M G GL^-^ (C) such that {Ay^ \ By^) is 
row- equivalent to {uo (M — I)\M ^ I), then rank {Ay^ \ By^) = dy^ . 



Ay^ and By^ serve only to represent the system of equations Ay^ • / 
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Remark. We recall from section [2] that this condition holds for all k when F's 
Laplacian is self- adjoint. Therefore, in this case ^/r is exact, as stated in theorem 

m 

Proof. Denote v = Vk, v = Vk, and recall that Ey = {gie^'}^^i is the set of edges 
entering v. Assume, by reordering if necessary, that Vi = /i^ for 1 < i < m, i.e., that 
{gie^'}^i are representatives for the G^-orbits in Ey. Denote = giC^' = gie^^ 
(where 1 < i < m), and note that G^i is conjugate to Gg^i . The action of Gy on 
Ey gives rise to a representation C [Ey] of G^, and the G^-set isomorphism Ey = 
Ili^i Gy ■ e'^ = iXiLi ^^//c-i translates to an isomorphism of G^-representations: 

m m 

C[£.]-0CfV/Gj-0Indg;lG,.. 

i=l i=l 

Here Iq denotes the trivial representation of a group G, but we shall also use it to 
denote its character. We now see that 



(4.6) {xc[e,]^Xr) 



G^ 



= ^(lG^_^,Xi?)^^^ =^dimi?^e~^ =^dimi?^^-^ =Y,d^,^=dv• 

i=l i=l i=l i=l 

We return to the matrices {Ay \By) ^ Mndx2d^ (C) and {Ay \By) e Mnx2n (C). For 
/ G (r), the action of Gy on Ey induces a permutation action of Gy on the entries 

of/ 



(f{v)] , and exactly the same action is induced on the entries of 



. Thus, the space C^^ of possible values and derivatives at v has naturally the 

structure of the G^-representation C [^^]0C [Ey]. Furthermore, as by assumption G 
preserves the boundary conditions, ker (^4^ \ By) C C^"^ is a sub-G^ -representation 
of C^^ = C [Ey] C [Ey]. We observe that the encoding and decoding processes 

are "rigid", in the sense that for x G [0,/e*] it suffices to know \ fr\ {x)\ 

I le* J reR 

and ^/| 



to 



determine < ^/ (x) > , and vi 
I ^5 J j=i..d. 



vice versa. Likewise 



ise, \fr\ I 



determine one another, and the same goes for the corresponding derivatives. This 
means that in the commutative diagram 



HomcG(i^,H(F)) 



HomcG, (i^,ker {Ay\By)) 



I 



(*/L'(*/)'L) 



■ker {Ay\B^) 



H{r/R) 

the map ip, which is this "local" encoding, is in fact an isomorphism. This gives us 
(4.7) null {Ay \ By) = {xR,Xker{A^ \ b,)) g~ ' 

so that by ([Q]) 

rank {Ay \ By) = {2xc[e^] - Xker(A, | b^)^Xr)g^ • 
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We therefore have 

rank {Ay \ By) = dy ^ {xc[e^] - Xker(A^ | b^), Xi?)^- = ^ ' 

and the last equahty holds for all representations of G if and only if Ind^.C [Ey] = 
Ind§_ ker {Ay \By). In particular, this happens if C [Ey] and ker {Ay \ By) are iso- 
morphic G^-representations, which we now show to follow from our assumptions. 

Observe that ^ : C [Ey] ®C[Ey] ^ C [Ey] , defined by ^ (a, b) = ua — bis di homo- 
morphism of G^-representations, and recall that ker {Ay \ By) is naturally embedded 
in C [Ey] C [Ey]. When restricting ^ to ker {Ay \By) we obtain the desired iso- 
morphism onto C [Ey], since dim ker {Ay \ By) = nuh {uj {M — I) \ {M + /)) = dy = 
dimC [Ey], and 

, 1 \ ( u;{M-I)a^{M^I)b = \ , ^ 

(a,6)GkerK ]^{ In >=>{a,k) = 0. 

□ 

4.3. Remarks. 

4.3.1. It seems of some interest to point out that the encoding process we have 
described has actually nothing to do with eigenfunctions of the Laplacian. The 
assumption that G (F) was not used during the construction of the quotient, 
and as a result, if no dummy vertices are introduced at the beginning of the con- 
struction, then we actually have 

(4.8) ^ : HomcG [R. C°^{^)) ^ C^i'^/R) • 

If dummy vertices are added, and T' is the graph obtained from T by their introduc- 
tion, we obtain only C°^{^/r) = HomcG (i^,G°°(r)), and unfortunately G°°(r) 7^ 
C (r), as was remarked in section [2l We have introduced dummy vertices in order 
to avoid loops and parallel edges, and also to ensure that a subset of the edges can 
be taken as a fundamental domain for the graph. Of these causes, only the last is 
unavoidable; one can still carry out the construction (encumbering somewhat the 
notations), even with loops or parallel edges, both in the original graph and in the 
quotient. The only case in which the construction fails altogether, and a dummy 
point must be introduced, is when a group element inverts the direction of an edge, 
so that a fundamental domain must include only half of the edge. Thus, in order 
for it to be possible to construct by our method a "smooth quotient", in the sense 
of (|4.8p . fix^r must be a subgraph of T for every g e G. We shall return to these 
observations in section 16.11 

4.3.2. If G acts on F and is a representation of H < G, we can consider the 
composition of isomorphisms 

n (r/i?) ^ BomcH {R. H {T)) i HomcG (indgi?, H (F)) ^ H (r/indgi?) 

where ^ and ^' are the isomorphisms defined during the constructions of ^/r and 
r/indgi?, respectively, and is the Frobenius isomorphisnj^. We obtain what is 
known as a transplantation (see \^Ax i25j) between ^/r and r/indgi?, an operator 

^Taking the induction to be the scalar extension Ind^i^ = CG^cH R, ^ is defined for decom- 
posable tensors by \Tf\ = 9 • fr (where / G Homci^ (R, H (F)), g ^ G, r ^ R), and extends 
linearly to all of the tensor product. 



14 



ORI PARZANCHEVSKI^ AND RAM BAND^ 



which construct functions on one graph as hnear combinations of segments of func- 
tions on the second graph. This is developed in more details in [14J. 

4.3.3. It is natural to ask, for a graph T whose Laplacian is self-adjoint, whether 
the Laplacian on ^/r is self-adjoint. This turns out to depend on both the action of 
G and the choices of bases in the construction, and it is addressed for some cases 
in [g. 

4.3.4. Another natural question is the following: for a quantum graph T acted 
upon by G, when does an irreducible representation 5 of G appear in 

It is known that every quantum graph with edges whose Laplacian is self- adjoint 
has a nonempty spectrum (see for example [23]). Therefore, if ^/s^s Laplacian is 
self-adjoint then S appears in H (F) iff ^/s has edges, and by the construction 
method this happens iff for at least one edge e in F the representation Res^^^ has 
a nonempty trivial component, i.e., (xs',!)^'^ ^ 0. In particular, if F's Laplacian 
is self adjoint, and G acts freely on F, then a self- adjoint quotient can always be 
obtained [14J, and each stabilizer has only the trivial irreducible representation. 
Thus, every irreducible representation of G appears in H (F). 



We now demonstrate several applications of the theory presented above which 
yield isospectral graphs. All the examples below are direct consequences of the 
theorem or the corollary presented in section [3l 



Figure 5.1. A graph that obeys the dihedral symmetry of the 
square. The lengths of some edges are marked. 



Let F be the graph given in figure [5Tl The lengths of the edges are determined by 
the parameters a, 6, c and it has Neumann boundary conditions at all vertices. G = 
D4, the dihedral group of the square, is a symmetry group of F. Denote by r the 
reflection of F along the horizontal axis and by a the rotation of F counterclockwise 
by ^/2. Then we can describe G and some of its subgroups Hi^ Hs < G by: 



5. Examples of isospectral quantum graphs 




2b 



G = {e, <j, <j^, cr^, r, r<j, ra^, ra^} 

Hi = {e, r, rcr^, a^} 

H2 = {e, rcr, rcr^, a^} 

Hs = {e, Cr, C^^ Cr^} 



This question, in the context of compact Lie groups acting on Riemannian manifolds, is 
addressed in [22j . 
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Consider the following one dimensional representations of i^i, H2 and respec- 
tively: 

(5.1) Ri : {e^(l), r ^ (-1) , ra^ ^ (1) , ^ (-1) } 

(5.2) R2 : {e^(l), ra ^ (1) , ra^ ^ (-1) , ^ (-1) } 

(5.3) i?3 : {eK^(l), a (z) , (-1) , (-z) } 

These representations fulfill the condition in corollary [H Ind^^i^i = Ind^^^s — 
Ind^gi^s and thus we obtain that ^/Ri^ ^/R2 and ^/rs are isospectral (figure [5^2]) . 




(a) (c) 

Figure 5.2. The three isospectral graphs ^/ri,^/r2,^/r3. Neu- 
mann boundary conditions are assumed if nothing else is specified. 
D stands for Dirichlet boundary conditions and N for Neumann. 

We now explain the process of building the graph ^/Ri. First we give an intu- 
ition which suffices to obtain the quotient in this case, and afterwards we strictly 
implement the method that is described in section [421 Going back to (|3.4p . we 
observe that the r.h.s. of it is 

(5.4) Home//, {Ri,n{V)) ^ {H{T)f' , 

where (7Y(F))^^ is the i?i-isotypic component of HiT) (considered as a CHi- 
module); the isomorphism is due to the fact that Ri is one-dimensional, hence 
irreducible. Let us study the properties of / G (7Y (r))^\ We know (see (|5.ip ) that 
rf = — /, which means that / is an anti-symmetric function with respect to the 
horizontal reffection. We deduce that / vanishes on the fixed points of r (marked 
with diamonds in figure [531 a)). 



1 o 

(a) ^ (b) 

Figure 5.3. (a) The information we have on f e (H(F))^\ Di- 
amonds mark the vertices on which the function vanishes and 
squares the vertices with zero derivative, (b) The quotient graph 
which encodes this information. D stands for Dirichlet bound- 
ary conditions and N for Neumann. 
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In a similar manner, we see that / is symmetric with respect to the vertical 
reflection since ra'^f = /, and therefore the derivative of / must vanish at the 
corresponding points (the squares in flgure [Of a)). Furthermore, it is enough to 
know the values of / restricted to the flrst quadrant (the bold subgraph in flgure 
I5.3r a)) in order to deduce / on the whole graph, using the known action of the 
reflections, which follows from f e {H (F))^^: 

(5.5) Tf = -f,Ta'f = f. 

Our encoding is now complete and the quotient ^/Ri is the subgraph which lies in 
the flrst quadrant, with the boundary conditions of Dirichlet and Neumann in the 
appropriate locations as was found for / (flgure [53Kb)). The encoding is described 
by the map ^ : HomcG (^i,^(r)) H{^/Ri) which is just the restriction map 
of functions in HomcG {Ri^'H (F)) = H (F)^^ to the mentioned subgraph. An im- 
portant observation is that given / G HiJ^/Ri) it is possible to construct a unique 
function / G (7Y(F))^^ (using (|5.5p ). whose restriction to the flrst quadrant sub- 
graph is /. It follows that ^ is invertible and thus is an isomorphism. This ends 
the intuitive approach and we now proceed to the rigorous derivation. 




(a) 

Figure 5.4. (a) The graph F with the representatives of 
^///i, '^/hi marked in bold, (b) The resulting quotient ^/Ri. 

First, we add "dummy" vertices to the graph F so that no vertex is carried 
by the action of Hi to one of its neighbors, and choose representatives je*}.^-^ 
for the orbits ^///i, and {vk}\^i for the orbits ^/hi. These representatives are 
marked in flgure I5.4r a) by bold lines and points. The dummy vertices amongst 
the representatives are i)i, i)2, ^5, ^6- ^i is one dimensional, and di = 1 for all i 
since the stabilizers of all edges are trivial. Therefore, the quotient graph is formed 
by taking one copy of each of the representative edges (flgure [5^ b)). Now, let us 
determine the boundary conditions using (|4.4p , (|4.5p . For all vertices we have d = 1 
and therefore Ay^ ^ Id = ^ = . Consider the vertex Vk = vs for 

which 

n = 3, m = 3, dy^^ = 3 
gi = e, i^i = /ii = 1, g2 = e, J^2 = 1^2 = 2, ^3 = e, = /j.3 = 3 
(5 = h 6 = (6' Id) =& = h 
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Plugging all this into (|4.4p , (|4.5p and using the boundary conditions on v^, which are 
given by Ay^^ ~ ( q o V ) ' ~ ( i i i ) ^^^^^ Neumann boundary conditions for 



'^3 as well: Ay 



A~ 



By^ = By^. Exactly the same treatment can be done for 



the vertex ^4 and the same boundary conditions are obtained. The case is different 
for the vertex ^'5: 



n = 2, 
gi = e, vi= 111 



1 



m = 1, 

= 4, 92 =T, = Ml = 4 

e = (e'0/,) = 6' = (l) 



The boundary conditions on v^, are of Neumann type as well: Ay 
{11). This time we obtain 



^ 1 -1 ^ 
.0 . 



By 



By 



= (??) 



'i)-(}) = 



(8). 



Ay^ and By^ are then reduced to square one dimensional matrices as expected, by 
removing the second row in both of them. We remain with Ay^ = (2), By^ = (0) 
which means Dirichlet boundary conditions on the vertex ^5. The same boundary 
conditions are obtained for vq. Similar derivation for vertices vi, V2 gives Neumann 
boundary conditions for each one of them. The rigorous construction thus gives us 
the same quotient graph that was obtained by the intuitive method (figures I5.2r a) , 
iLlb)). 

The quotient ^/r2 can be constructed in a similar manner, and is shown in fi.g- 
ure I5.2r b). We proceed to demonstrate the construction method for the quotient 
V^sLJ. We first add the corners of the square as dummy vertices to F {vi in figure 
I5.5r a) is one of them). We are not obliged to do so, but it yields a quotient with 
simpler boundary conditions. The choice of representatives for the edges and the 
vertices is shown in figure [531f a) and the resulting quotient in figure [531 b) . 
The vertices V2 and ^'3 have Neumann boundary conditions exactly as their pre- 



ae 




(a) 




By 



(b) 



Figure 5.5. (a) The graph F with the representatives of 
^/hs, ^/h3 marked in bold, (b) The resulting quotient ^/r3- '^2, '^3 
possess Neumann boundary conditions. 



-'^-'^This result was obtained with G. Ben-Shach. 
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decessors, V2 and v^,. For vi we obtain more interesting boundary conditions: 

^vi — \o / ' -'^'^1 ~ V 1 1 / 

n = 2, m = 2, dy^ = 2, 
gi = e, i^i = /ii = 1, 92= CF, = /i2 = 4 

which gives 

(5.6) A,,={li), = (?oj. 

Non-formahy speaking, the vertex vi "apphes a factor of z" to the functions that 
cross it. The resulting graph is the one that was shown in figure lOlf c) . 

In order to exhaust this example, we observe that Ind^^i^i = Ind^2i?2 — 
Ind^gi?3 is the two-dimensional irreducible representation of 1^4, which we denote 
by R. By theorem [3l the isospectral family of the three graphs given in figure 15.21 
can be extended by adding any graph which is ^/r. We therefore construct now such 
a graph. Let us use the intuitive approach first. Recall that (|5.4p was the key for 
the intuitive construction of ^/Ri. Analogously to (|5.4p . we make the observation 
that encoding Homcc (F)), the r.h.s. of (|3.4p . is similar in nature to encoding 

(7Y(F))^, the i?-isotypic component of H(F), as due to the simplicity of as a 
CG-module the two are isomorphic. This can be understood as follows: making a 
choice of a basis 62} for i?, and given a function / G HomcG {R^ ^ (F)), we have 
that G (H (F))^ and furthermore j/bi, fh2 1 spans over C a CG-module iso- 

morphic to R. In order to exhibit the general behavior we avoid sparse matrices, 
and pick a basis {61, 62} for which the matrix representation of R is 

It is enough to consider only the matrices of these two elements for the construction 
of the quotient. 




(a) 

Figure 5.6. (a) Two copies of the graph F with the representa- 
tives of ^/i^4, marked in bold. These two copies are merely 
a visualization of the "basis functions" fbi^fb2 on F. (b) The first 
stage in the formation of ^/n is the gluing of both copies in the 
vertex V4, with the boundary conditions given in (|5.12p . (|5.13p . 
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Examine the properties of /b^, that fohow from the above matrix represen- 
tation (figure [5?6lf a)). Since / G HomcG {R,H{T)) we have ra^fb^ = /^^^s^-i^^ = 
fra^bi^ ^iid thus the first column of the matrix representing ra^ tehs us that 

ra'k = ^/2fb,-y2h, 



(5.8) 
(5.9) 



and enables us to relate the values and the derivatives of Z^^, fb^ on the vertex V4. 
Since V4 is a fixed point under the action of ra^ and there are Neumann boundary 
conditions on it, we have that 

(5.10) (ra'k) {V4) = k (va) 

\ / 63 63 

(5-11) (rcT^fiMA^i) = -4L(«4). 

\ / Ie3 1 63 

Evaluating (|5.8p on V4 and combining this with (|5.1Qp gives 

(5.12) (l-V3/2)^ (^4)+V2/6. (^4)=0. 

63 63 

Similarly, from (|5.9p and (|5.1ip we obtain 

(5.13) (-l-v/3/2)4 (f}4) + V2/6, (^)4)=0. 

63 63 

We may therefore think of two copies of the graphs. Each of the basis functions 
/fo^, resides on one of the copies, and the relations between the values and the 
derivatives of the functions allow us to take a subgraph out of each copy (marked in 
bold in figure [5^61 a)) and glue both of them together with the appropriate boundary 
conditions. The first stage in this gluing process, visualized in figure ISTBlf b). is to 
identify the vertex in the two copies and turn it into the vertex V4 of the quotient 
with the boundary conditions that were derived in (|5.12p . (|5.13p : 



(5.14) 



( I-V3/2 1/2 

V 



f 0\ 

y -1-V3/2 1/2 J • 



After treating similarly vertices vi^V2 we get the quotient ^/r (figure [5?7|) whose 
remaining boundary conditions are given by: 

(5.15) A^, = = ( t 7^ ) , = B,, = ( 4, ) . 




Figure 5.7. The quotient graph ^/r which is isospectral to the 
graphs in figure [5^21 The boundary conditions are as described in 
([5:T4I) . ([5T51) . 
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We now use the rigorous approach for the same quotient, ^/r. The represen- 
tatives of the orbits ^/g are {e*}^^-^ and the representatives of ^/g are {vk}\^i 
(figure [SSJ a)). This time the representation is not one dimensional {d = 2) so 
there are additional details to consider. First, note that we have two copies of each 
representative of ^/g in the quotient and both of the copies survive since all edges 
have trivial stabilizers (figure ISTSlf b)). This last observation ensures that we can 
take = B for all i (i.e., the same basis for all edges). We again take B to be the 
basis for which the matrix representation of R is (|5.7p . 





(b) 



Figure 5.8. (a) The graph F with the representatives of 
^/d4, ^/d4 marked in bold, (b) The resulting quotient ^/r. 



We treat the boundary conditions at the vertices one by one: 
• ^4 has the following data: 



m 



1, 



6 



91 



e, z/i 



^1 = 

(1)0/2 



3, 



92 



■ TCF 



^2 = [^\ 



1 

1 

1 
1 



, Br^^ are the regular Neumann matrices and we therefore obtain 



Bv 





Noting that both Ay^ and By^ are of rank one, we see that they express the 



same boundary conditions as given in (|5.14p . 

vi obviously has the same boundary conditions as V2- We examine ^i: 



2, 



1, d„ 



91 



Vi = 111 



1, 

1 
- f 1 

1 
1 



92 



^2 = Ml 



Again, Ay^ and By^ are the regular Neumann matrices and we get: 
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which are matrices of rank one and again we may reduce these matrices 
into two dimensional ones which are exactly those given in (|5.15p . 
• The case of is a bit more interesting: 



gi = e, z/i = /ii = 1, 



n = 3, m = 3, dy^ = 6, 
e = (O' h) = h^l2=h 



As Ay^ and By^ 



are Neumann matrices, we have 



Ay 



By 




and we see that the above boundary conditions separate the edges into two 
sets, {ex,e2,ef} and {62,62,62}, each dominated by a regular Neumann 
condition. This enables us to split the vertex ^3 into two distinct vertices of 
degree 3, each connected to a different set of edges and possessing Neumann 
boundary conditions. We remark that this would happen for any choice of 
basis for as here gi = 92 = 93 = 

Note that the resulting quotient is the same as was obtained previously (figure [5?7|) . 

Finally, we repeat the construction for an arbitrary choice of basis which yields an 
orthogonal matrix representation for R. We can parametrize such a representation 
in the following way: 



ra 



cos U — sm U 
—2 cos ^ sin ^ 

2 cos sin 

.2, 



2 cos sin 

2 n , .4^2 , 



COS" 

.2 



sm 

2, 



cos 



sm 



cos u — sm 
—2 cos ^ sin ^ 



For example, the basis we chose in (|5.7p is obtained hj 6 = ^/s. As remarked, 
V3 always splits into two vertices with Neumann conditions, so that figure [5?7l can 
describe the quotient with respect to any basis. For the parametrization above, we 
obtain the following boundary conditions: 



A —A 



f 2sin^ 9 sin 26* 
sin 20 2-2sin^0 

V , 

^ 


2-2 sin^ 6 - sin 29 
- sin 29 2sin^ 9 , 



= 



l-sin26' 2sin2 6'-l 
2sin2 6'-l l+sin26' 





l+sin26' l-2sin2 6' 
l-2sin2 6' l-sin26' 



All of these matrices are of rank one, and can therefore be reduced to square ones 
by deleting the appropriate row£l. We thus get a continuous family of isospectral 
graphs. Examine two members of this family: = and = ^^4. The boundary 



■"^ ^However, there is no a priori reduction which is vahd for all ( 
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conditions for the case = are: 

A —A — (0 2\ A — (1 -1\ 

B — B — f20\ o — fOO\ 

J^vx — -Dv2 — 1 y -^V4 — \1 1 ) 

When applying this to figure 15. 7[ we notice that the vertices vi , V2 do not stay 
vertices of degree two, but rather, each of them sphts into two vertices of degree one, 
one with Dirichlet boundary condition, and the other with Neumann. The vertex 
V4, however, stays connected and obtains Neumann boundary conditions. Observe 
that the resulting quotient is the one that we have already obtained as ^/ri (figure 
I5.2r a)). In a similar manner, the quotient ^/r2 (figure lOlf b)) is obtained from the 
choice 6 = ^^/a. We conclude by pointing out that the graph described in figure 
15.71 is a good prototype for the mentioned isospectral family, yet it might also be 
misleading, since there are members of the family whose boundary conditions tear 
apart the edges connected to some of the vertices and thus change the connectivity 
of the graph. One should also pay attention to the fact that we have treated only 
orthogonal representations of D4. These are not the most general ones, and we 
may extend the isospectral family presented above by considering the broader case 
of all matrix representations of R. In particular, the quotient ^/Rs (figure lOlf c)) is 
obtained from the unitary representation 

-0' ^^(-1 "0') }• 

6. Isospectral manifolds and stratifolds 

If r is a Riemannian manifold equipped with an action of a finite group G, 
then C°°(r) is again a module over CG [x] , with x acting as the Laplace-Beltrami 
operator A. If however F has a boundary, at which differential boundary conditions 
are imposed on C°°(r), then in general it is no longer closed under A. In order to 
treat this case as well, we limit our attention to the subspace of C°°{T) spanned by 
A's eigenfunctions, which we again denote by H (F). Assuming that the boundary 
conditions are linear, H (F) is closed under A and is therefore a CG [x]-module as 
before. Section [3] is naturally generalized to these settings: 

• For a representation R of G, we define a r/i?-manifold to be a Riemannian 
manifold (possibly with boundary, at which homogeneous conditions are 
imposed) F', such that there is an isomorphism 

(6.1) H(r) = HomcG(i^,H(F)) 

intertwining the Laplace-Beltrami operator. 

• For a representation RoiH <G^^/r and r/indgi? are isospectral; therefore, 
r/cG is isospectral to F, and for representations Ri^R2 of i^i,i^2 < G 
satisfying lnd%^Ri = Ind^^^s, ^/Ri and ^/r2 are isospectral. 

The main advantage of quantum graphs for our purposes is that under fairly mod- 
erate assumptions (e.g., self-adjoint Laplacian or a free action) one can build a 
quotient for every representation, as is demonstrated in section [H 

Graphs are one-dimensional manifolds with singularities (at the vertices), and 
it is these singularities that we exploit, by endowing them with the appropriate 
boundary conditions, to encapsulate the restrictions arising from a choice of a rep- 
resentation. In higher dimensions, manifolds with a boundary, carrying Neumann, 
Dirichlet, or a more complicated boundary condition, are a generalization of this 
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idea, and one goal of this section is to demonstrate that some known isospectral 
examples of such objects can be understood using our theory. That is, we show that 
for some known isospectral pairs the manifold and boundary conditions are such 
that the objects are quotients (in the sense of definition [1]) of a common manifold by 
two representations with isomorphic inductions in some supergroup of symmetries. 

It turns out, however, that in order to form a quotient by a general representation 
we need more singularities than just boundaries (at least via our construction). A 
graph is a one dimensional manifold when all of its vertices are of degree two, and a 
manifold with boundary when all vertices are of degree at most two. Unfortunately, 
even if a graph has one of these properties, its quotient by a multidimensional 
representation (as constructed in section [4]) need not have either, since the degrees 
of the vertices are multiplied, in general, by the dimension of the representation. 

Carrying over the construction method of section [4] to general Riemannian mani- 
folds (e.g., by replacing graphs with higher dimensional simplicial structures) yields 
objects we might call "quantum-stratifolds". In general, these consist of several 
Riemannian manifolds of the same dimension "glued" along their boundaries by 
homogeneous boundary conditions (so in dimension one, we obtain the notion of 
quantum graphs). When a boundary condition involves the boundaries of more 
than two manifolds, the result is no longer a manifold, but rather a stratifold. 
Even though this is in general the case, by choosing an appropriate action, repre- 
sentation and bases, it is possible to obtain manifolds even when taking a quotient 
by a multidimensional representation. 

6.1. Isospectral drums. In [26t [27]. Jakobson et al., and Levitin et al., respec- 
tively, obtain several examples of isospectral domains with mixed Dirichlet-Neumann 
boundary conditions, all of which can be interpreted as quotients with respect to 
representations sharing a common induction. As a basic demonstration of the gen- 
eralization of our theory to higher dimensions, we reconstruct an isospectral pair 
consisting of a square and a triangle with mixed boundary conditions (figure 1 in 



Figure 6.1. The two isospectral domains presented in [27], ob- 
tained as quotients of the square S (figure 16. 2p by the representa- 
tions in (|5.ip . (|5.2p . Solid lines indicate Dirichlet boundary condi- 
tions and dotted ones Neumann. 

This example rests upon our acquaintance D4, so that we can reuse the defini- 
tions and results of section [H In place of the graph in figure \bA\ we now consider 
the full square S, with Dirichlet boundary conditions, and with G = D4 acting as 
one would expect (figure [6^2]) . 

The domains in figure 16.11 are quotients of the square S (figure 16. 2p by the rep- 
resentations Ri and R2 of i^i,i^2 < G, which are defined in (|5.ip . (|5.2p . Since 

G G 

Indji^^i?! = lnd[f^R2, the two domains are isospectral. 



[22], O here). 
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rcr 



T(J 



Figure 6.2. The square S, and the axes of reflection elements in D/^. 



We demonstrate the construction of ^/ri. Recahing that Homc//i (i?i,H(S)) = 
(7Y(S))^\ we have again that Hi^/Ri) should encode the -isotypic component 
of l-t{S). T, the flrst quadrant of S (flgure lOlf a)). is a fundamental domain for 
the action of Hi , so that given f e 7i (T) it is possible to construct at most 
one function in (7Y(S))^^ whose restriction to T is /. Thus, the restriction map 
^ : (7Y(S)) ^ (T) is injective. In order for it to be surjective, we must 

impose suitable boundary conditions on T. From (|5.ip we obtain information on 
/ G (7Y(S))^\ Since such / is anti-symmetric with respect to the action of r, 
it must vanish at the horizontal axis of reflection, and therefore every / G im^ 
vanishes at the lower edge of T. Similarly, every / G (7Y(S))^^ is symmetric with 
respect to rcr^, so that its normal derivative at the vertical axis of reflection is zero, 
and thus all functions in im ^ have vanishing normal derivatives at the left edge of 
T. This information, summarized in flgure iGTSlf a). suggests the domain presented 
in flgure lOlf b) as the quotient ^/ri: a square identical to T, three of whose edges 
have Dirichlet boundary condition and one Neumann. 



. rcr^ 1 



T 



r —1 



(a) 



(b) 



(c) 



Figure 6.3. (a) The fundamental domain T for ^/hi; every / G 
{H (S))^^ vanishes along the dashed line and has zero normal de- 
rivative at the dotted line, (b) The quotient planar domain ^/ri 
which encodes this information. The solid lines represent Dirich- 
let boundary conditions and the dotted one Neumann, (c) The 
square of which (b) is a smooth quotient; along the dashed lines 
functions need only be once continuously different iable. 



Once these boundary conditions are imposed on ^ is indeed onto: for 

/ G {^/Ri) which obeys them, we deflne a function / on S by / = /, t/ = — /, 
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rcr^/ = /, (T^/ = — /. While / is well defined on the vertical ra^-axis even if / 
does not obey any boundary conditions, it is the requisition that / vanish on the 
lower edge of T that guarantees that / is well defined on the horizontal r-axis. In 
a similar manner, while at the r-axis the two one-sided normal derivatives of / 
agree a priori, it is the Neumann condition at the left edge of T which ensures this 
at the rcr^-axis. The boundary conditions thus assure that / is well defined and 
continuously different iable, and being piecewise smooth and a sum of Laplacian 
eigenfunctions, it is smooth, and therefore in TL (S), so that / = ^/ G im^. As ^ 
and its inverse are obviously C [x]-linear, we have established Homc//i {Ri^H (§)) — 
{n ^ n as the definition of a §/i?i-domain in 1^ calls for. 

Analogously, from the properties of / G (7Y (S))^^ we can deduce the correspond- 
ing quotient ^/i?2. This process is summarized in the two parts of figure lOl 

_ ra 1 



/ 

/ 

/ 

''/■. 

/ 

/ 

/ 

/- 

/ 

(a) 

Figure 6.4. (a) The information we have on / G (H(§))^': it 
vanishes along the dashed line and has zero normal derivative at 
the dotted line, (b) The quotient planar domain ^jRi which encodes 
this information. The solid lines represent Dirichlet boundary con- 
ditions and the doted one Neumann. 

We return once more to the question of smooth quotients, which was raised in 
section [4.3.11 It turns out that even with the boundary conditions we have imposed 
on ^ is not surjective as a function from ^ to Consider 

for example the smooth function x^'ip' ^ when regarding the lower-left corner of T 
as the origi It is not the restriction of any function in (C°°(S)) \ It is the 

restriction of a function in (C^(S))^\ namely, x^y • \y\. Once again we encounter 
the problem of modeling smooth structures by differential boundary conditions of 
degree one, which was addressed in section [2l In fact, ^ establishes an isomorphism 
^ (C°°(§0)^\ where is § after the addition of two "Neumann lines of 
degree two" at the axes of refiection corresponding to r and ra^ (figure lOlf c)): it 
is a square divided into four, such that a smooth function on is by definition a 
function which is smooth on each inner (closed) square, and satisfies the Neumann 
condition where two squares meet, or equivalently, is at the r and tg^ axes. 
Had we proceeded by the rigorous method presented in section [421 we would have 
added these Neumann lines before the construction, in order to allow a fundamental 



rcr 



rbi 



It is not really a function on T, since it does not vanish at the right and upper edges. This 
can be rectified by taking x^ip' (x^ — l) (y^ — l) instead, but it would clutter the argument. 
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domain which is a subcomple S A gain, this could be justified by the preservation 
of spectral properties: we have, as for graphs, that H{E>) = 7Y(S'), so that a ^R- 
quotient (now in the sense of definition [1]) is the same thing as a ^'/i?- quotient. 

We remark that the various constructions demonstrated in section \5\ can be ap- 
plied analogously to S, enriching the isospectral pair in figure 16.11 For example, 
^/r3 would be an orbifold with a line that applies a factor of i to functions crossing 
it. The other isospectral families in [26^ J7J can be obtained from various repre- 
sentations of the general dihedral groups Dn, and of the product D4 x D4. The 
interested reader will find some of these constructions in fl4]. 

6.2. The Gordon- Webb- Wolpert drums. In a similar fashion, we can apply our 
method to the Gordon- Webb- Wolpert construction [31 [4], obtaining their isospectral 
planar domains with new boundary conditions. We follow the exposition of Buser 
et al. [5], who obtain the mentioned drums as follows: they consider Go, a group 
of motions of the hyperbolic plane H (*444 in Conway's orbifold notation), and an 
epimorphism n : Go ^ G = PSL3(2). In G they exhibit two subgroups A and 
5, each isomorphic to ^4, that satisfy the Sunada condition [Ij with respect to 
G. The quotients of H by 7r~^{A) and 7r~^{B) are isometric domains. Both are 
composed of seven copies of a hyperbolic triangle (which is a fundamental domain 
for the action of Go), assembled in different configurations (which are determined 
by the coset structure of the pre-images). Finally, by replacing the fundamental 
hyperbolic triangle with a suitable Euclidean one, the non-isometric isospectral 
drums of Gordon et al. are obtained. 

An elegant formulation of the Sunada condition for Hi and H2 in G is that the 
inductions of the trivial representations Ihi and 17:^2 ^ isomorphic, i.e. 

(6.2) Indg^lj/, ^ Indg^l^^, . 

In fact, the connection between A and B is stronger than this (refiecting a line- 
point duality in the Fano plane): it turns out that for every representation R 
of ^4, Ind^R ^ lnd%R. For each such R, we can thus construct an isospectral 
pair by taking the quotient of EI by the pullbacks of R to 7t~^{A) and 7t~^{B). 
Taking R = Is^ will produce once again the planar drums of Gordon et al. In 
fact, we shall see in section [631 that taking quotient (in our sense) by the trivial 
representation of a group is equal to taking quotient (in the classical sense) by 
the group. Taking R to be the sign representation of ^4, and again replacing the 
fundamental hyperbolic triangles with Euclidean ones, we obtain the same drums 
but with different boundary conditions (figure [63]). 

We conclude this example by pointing out that in [5j a wide variety of isospectral 
pairs is presented, using various symmetry groups of H. All these examples can 
be exploited to construct other isospectral pairs, as isomorphic inductions may 
be found either from Sunada triples or by taking appropriate sums of irreducible 
representations. 



Each pair of opposite sides in S is swapped by some element in Hi, and as in dimension 
one we did not allow a vertex to be moved to a neighbor by a group element, neither should be 
lines allowed to, in dimension two. In a more general view, for FixrS and Fix^^2S to be each a 
subcomplex of S, S must be subdivided (by Neumann lines) into S^ 
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(a) (b) 

Figure 6.5. The isospectral drums of Gordon et al. with new 
boundary conditions. 



6.3. The Sunada method. We recah the classical theorem of Sunada [T]: 

If G acts freely on a Riemannian manifold F, and i^i,i^2 ^ G 
satisfy (|6.2p . then ^/hi and ^/h2 are isospectral manifolds. 

Sunada's theorem follows from the definition and corollary at the beginning of the 
current section, once we show that for a finite group G acting freely on a manifold 
F, the quotient manifold t/g is a r/i ^-manifold, that is, 

(6.3) n (r/G) ^ HomcG H (F)) . 

This follows from the observation that Homcc (ic, C'°°(F)) corresponds natu- 
rally to C°^{T)^^ = C°°(F)^, the trivial component of C°°(F), and this is the space 
of functions on F which are stable under all elements of G. But these are exactly 
the functions which factor through t/g, hence C°°(r/G) = HomcG (Ig, ^°°(r)), and 
in particular (|6.3p follows. 

Remark. We can view the preceding argument as yet another proof for Sunada's 
theorem, but this would be presumptuous. In fact, Pesce ^28j uses Frobenius Reci- 
procity in exactly the same manner to reprove Sunada's theorem. A survey of 
different proofs for Sunada's theorem, among them Pesce's, can be found in [9]. 



7. Summary and open questions 

The main construction presented in this paper is that of objects denoted ^/r, 
where is a complex representation of a finite group acting on a geometric object 
F. For such F and R there can be, in general, many objects so denoted, and they 
are all isospectral to one another. Furthermore, these objects are defined so that 
whenever Homc//i {Ri^ _) — Homc//2 (^27_), where each Ri is a representation 
of a group Hi acting on F, there is also an isospectrality between ^/ri and ^/r2. 
The consequences of this are explored in section [3l and in particular we find two 
convenient means for the construction of isospectral objects: 
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Starting with a group G: take subgroups i^i,i^2 < G and corresponding 
representations Ri , R2 sharing a common induction in G For any object T on 
which G acts by symmetries, ^/ri and ^/r2 are isospectral. 

Starting with an object F: find a group G acting on F and construct ^/cG 
(by some choice of representatives and bases, as explained in section [42]) . Any 
quotient thus obtained is isospectral to F itself, by the analogue of proposition [2] 
for arbitrary geometric objects. 

It is natural to ask to what extent the various methods for obtaining isospectral 
objects overlap. For example, in section [5l three isospectral graphs (figure [5^2]) are 
obtained from representations with isomorphic inductions, but at the end of the 
same section it is demonstrated that all of them (together with others) could have 
also been obtained as ^/r for a single R (by different choices of bases). Can one 
expect that given a basis for i?, there is always a basis for lnd%R with respect to 
which ^/r and ^/ind%R are isometric? 

Even when limiting to the basic quotient construction, questions arise. For R 
and F as above, we have a family of isospectral objects ^/r^ varying as one moves 
between different choices of bases in the construction, as explained in section 14.21 
and demonstrated in the last part of section [5l This family has the topology of a 
manifold, being parametrized by the action of a general linear group on the space 
of possible bases. Surveying this continuum of quotient objects, one might ask 
where along it occur changes in the shape of the objects (in contrast with only 
boundary condition changes), in the number of connected components, etc. One 
can look for certain types of objects in this continuum, such as manifolds, billiards, 
objects with real boundary conditions, or ones with a self-adjoint Laplacian. Such 
questions seem to lead to a deeper research in differential and algebraic geometry, 
investigating the critical points at which changes occur or the algebraic varieties at 
which certain conditions are fulfilled. Except for the basic demonstration of these 
phenomena in section [5l we have not treated these questions. 

We list some more questions that seem interesting, and which we have not re- 
garded: 

• F is naturally a r/cG-graph. Does it occur by our construction? It seems 
that the answer is yes, by taking G as a basis for CG, but we have not 
shown this. 

• Given two isospectral objects, can it be decided algorithmically whether 
they are representation-quotients of a common object? 

• What are the necessary and sufficient conditions for the quotients con- 
structed in section [42] to be proper quotient graphs (in contrast with gen- 
eralized ones)? Exact quantum graphs? Graphs with a self-adjoint Lapla- 
cian? 

• Can the isomorphism (|3.4p be made natural, in a suitable category? This 
can be interpreted both as (contravariant) functoriality in R^ or as functo- 
riality in F, which would require a definition of quantum graph morphisms. 



^This resembles the Sunada condition, but is dramatically easier to achieve, since we are 
free to take any representations of the subgroups (instead of only the trivial ones). A systematic 
approach would be to take all irreducible characters of subgroups of G, induce them to G, find 
linear dependencies, and sum the corresponding representations accordingly. Also, any Hi and 
Ri are usable with H2 = G, by taking R2 = Ind^^ Ri . 
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• Can the theory presented in this paper be apphed to discrete graphs? To 
representations of Lie groups acting on Riemannian manifolds? 

• It is clear that n{Tl}r) =H (F) © H (F'), so that ar y r' = err + err' , and 
given bases for R and R\ their union is a basis for with respect to 
which ^/r®r' is isometric to ^/rU^/r\ Is there an operation on graphs, 
or general geometric object, which gives (F (g) F') = H{T) (8)7Y(F'), so 
that ar^v = crp • crp/? What about convolution: or^r' = or T*rcrr'? 

• A classical conjecture, originally aimed at Riemannian manifoldq^j [29]: for 

r 

G = AutF, and R = where {Si}l^^ are the irreducible represent a- 

i=l 

tions of G, is (Tp < 1? 
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